UNCLASSIFIED 


LINEAR  SAVES  ESTIMATORS  OF  THE  POTENCV  CURVE  IN 
BIOASSAV(U)  COLUMBIA  UNIV  NEW  VORK  L  KUO  1984 
AFOSR-RR -87-1882  4AF0SR-87-8872 

F/G  8/1 


SC 


LTNCLAS  S I  FI  ED 


2*  S£C-*'rv  CAASS.PiCATiCN  ai.*-C«i”y 


oe  c-*ss;»i  cation. oc^nc*  a  oing  sc*  e  out. 


3  3iST«ieuTic*.'A vAiwAaic:”v  *e»c*T 


Approved  for  public  release;  distribution 


irrl  *  . 


l.  MGNiTcflJNG  Organization  *6*0*7  num8£*:S: 

AFOSR TR-  87-  10  83 


64.  SAME  OS  »£S«;PM.NS  C.°Oa.niCa~iOn 

Columbia  University 


6*.  a£  C  *  £  35  f  C.  f>  Siczt  on a  C  jc# 

Columbia  University 
New  York,  NY  10027 


wNOiNG.'SPCNSCfliNG 

T:CN 


Sc  *GG*ESS  ;C.*>.  5  raw  an  a  Z.'?  ZoC*t 


A iNG  ORGANIZATION 


.cs  c:  Scientific  Research 


"b\c\  410 


Sc  CPPiCS  SYMSC.  9.  ®*»CCUS6MSNT  .NSTauMEsT  iC£NT,PiCAT;CN  NUM96  = 

ilf  «,,***«.  ftfosR  87-0072 


»“CCaAM  »*Cj£C"  !  -aSic  NCa*.  - s  1  ■ 

S.iWtNT  NC  |  %C  i  _  NC.  I  NO 


i  AS 


Beilin?  AFB  DC  ICI-Bl-ci-S 


y\  7i7„E  Insitmiii  iacu/'try  ?la*si f; nation 

Linear  Bayes  Estimators  of  the  Potency  Curve  in  Bioassay 


1J.  ®SaSONAu  AOT-CaiS. 


Lynn  Kuo 


*9.  ASSTAaCT  'Co*wnw«  on  nrt.#-»«  if  i«c(«i«a  ana  ic**a,'y  >>  nu m&ar. 


The  prior  distribution  on  the  class  of  potency  curves  in  quantal  bioassay  is  assumed 
to  be  Ferguson's  Dirichlet  distribution.  Given  the  integrated  squared  error  loss  and  the 
quantal  response  data,  we  derive  the  Bayes  estimator  in  a  linear  space  generated  by  the 
data.  Some  numerical  examples  and  asymptotic  results  are  also  given. 


DTIC 

lELECTE 

SEP8919W 


I 


rnnrr 


nam£  C*  *eS>£NS.Si*£  .NCiVtCUA^, 

3.  ’.vccdlruff 

DC  FCRM  1472,  22  APR 


22a  T£u5»v*CN6  nwm8E* 
ftnelua*  A 'tb  Co**  i 

(202)  767- 


"c  CP  P  tCS  SYMJOL 


E  3.  TiOY  C®  I  .AN  Y3  S  OaSCi.S”8. 


UNCL. 

SSS'-aiTY 


0 


t'JVJ 


!« 


m 


»%!«' 

i‘,V 

//«*: 

V,’ 


.  -»<• 

m 


AFOSR-TK-  87  -  10  83 


Linear  Bayes  Estimators  of  the  Potency 


Curve  in  Bioassay 


LYNN  KUO 


Division  of  Statistics,  University  of  California,  Davis 
Davis,  California  95616  U.S.A. 


SUMMARY 


Tne  prior  distribution  on  the  class  of  potency  curves  in  quantal  bioassay 


is  assumed  to  be  Ferguson's  Dirichlet  distribution.  Given  the  integrated 


squared  error  loss  and  the  quantal  response  data,  we  derive  the  Bayes 


estimator  in  a  linear  space  generated  by  the  data.  Some  numerical  examples 


and  asymptotic  results  are  also  given. 
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1.  INTRODUCTION 

The  Bayesian  nonparame cric  approach  Co  estimating  the  tolerance 

distribution  in  quantal  bioassay  has  received  some  attention.  The 

computational  difficulty  in  evaluating  these  Bayes  estimators  has  hindered 

their  applications.  This  paper  explores  the  linear  Bayes  approach  to  the 

bioassav  problem.  These  linear  3ayes  estimators  can  be  computed  easily  by 

using  statistical  software  which  has  the  capability  of  inverting  a  matrix. 

Let  us  state  the  quantal  bioassay  problem  as  follows:  The  experimenter 

intends  to  test  che  potency  of  a  stimulus  bv  giving  subjects  injections  of  the 

L  f,  u 

srimuius  at  different  levels;  namely,  he  chooses  L  dosage  levels,  t.,...,ty, 
i  *  b 

an c  treats  n . ,n  subjects  at  these  levels  respectively.  Each  subject 

X  t> 

possesses  a  fixed  tolerance  level.  If  a  stimulus  exceeds  a  subject's 
tolerance  level,  the  subject  responds  positively.  If  not,  there  is  no 

response.  Therefore  we  observe  the  number  of  positive  responses  at  each 

i  0  > n  u 

level.  These  numbers  are  denoted  by  The  potency  curve  F  is  the 

distribution  of  tolerance  levels;  i.e.  F  is  defined  by  the  probability  F(t)  of 
getting  a  positive  response  to  a  dosage  at  level  t  for  all  t.  The  objective 
of  this  article  is  to  make  inferences  about  the  potency  curve  F.  ■  • 

The  maximum  likelihood  estimator  of  F  subject  to  the  monotonicity 
restriction  may  be  found,  as  in  3a r low,  Bartholomew,  Sremner  and  Brunk  (1972), 
to  be  a  weightec  isotonic  regression.  If  we  have  some  prior  intormation  or 
F,  we  would  expect  to  do  better  bv  a  Bayesian  nonparametric  method.  Krart  and 
Van  Eeden  (1961)  derive  an  estimate  of  F  by  using  a  tailfree  process  prior. 

The  Bayes  estimator  of  F  using  Ferguson's  Dirichlet  process  prior  ( 1 9  -■  3 )  has 
seen  studied  by  Ramsey  (1972),  Ar.toniak  (1971),  Bhattacharya  (1981),  Disch 
(1981),  and  Amman n  (198-).  Ramsay  discusses  the  posterior  distribution  and 
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its  mode.  He  also  presents  interesting  examples  to  illustrate  the  desirable 
features  of  the  Bayes  estimators.  Bhattacharya  obtains  the  asymptotic 
posterior  distribution  by  means  of  a  Markov  Chain  formulation.  Ammann 
generalizes  the  Dirichlet  prior  to  a  class  of  priors  which  uses  processes  with 
independent  increments. 


Antoniak  (1974)  has  treated  this  problem  with  integrated  loss  function 
.2 


L(F,F)  =  j  (F-F)  dW  where  W  is  a  known  weight  function.  He  has  shown  that  the 
posterior  distribution  of  F  given  the  data  is  a  mixture  of  Dirichlet  process 
distributions,  and  he  has  derived  the  3ayes  estimator  of  F  for  two  dosage 
levels.  Unfortunately,  this  mixture  becomes  increasingly  intractable  when  the 
number  of  stimulus  levels  increases.  It  is  shown  by  Antoniak,  when  L  =  2, 
that  a  double  summation  (two  indices)  can  be  used  to  express  this  mixture; 
when  L  =  3,  six  indices  of  summation  are  needed.  In  general,  L(L-l)  indices 
are  needed  to  represent  this  mixture.  It  already  takes  a  substantial  amount 
of  bookkeeping  for  the  mixture  even  in  the  case  L  =  3 .  Therefore,  it  is 
important  to  search  for  methods  of  approximating  the  Bayes  rule. 

Disch  (1931)  proposes  two  methods  for  approximating  the  posterior 
distribution.  One  method  uses  a  single  Dirichlet  to  approximate  the  mixture 
of  Dirichlets.  The  parameters  of  the  single  Dirichlet  are  defined  to  be  a 


'unction  of  Z  k  and  Z  (r.  -k  ).  The  other  method  is  a  two-step  method. 
1  ..11 

Id 


T he  two  steps  involve  deleting  a  subset  of  observational  doses  and  using  the 
above  method  of  single  Dirichlet  approximation  again.  As  pointed  out  by 
jusch,  the  drawback  of  the  two  step  method  is  that  the  approximation  depends 


_  c  c  o  r  r\  r  A  r\  c  o  c  ■>  c  V  £>  *  ■"?  rr  *  of 


A  linear  Saves  estimator  is  proposed  here  to  approximate  the  Bayes 
estimator.  Given  the  loss  function  L(F,d)  =  /(F(t)-d(t))  dW(t),  we  restrict 
the  decision  rules  d(t)  in  the  linear  space  generated  by  k^,...,kL  and  1.  A 
linear  3ayes  rule  is  the  Bayes  rule  in  the  linear  space.  The  solution  is 
given  by  pointvise  minimization.  For  each  t,  we  find  constants 
depending  on  t,  to  minimize  E(F(t)  -  -  X^k,  -  ...  -  \TkL)“.  It  is  much 

easier  to  evaluate  the  linear  Bayes  rule  than  to  evaluate  the  Bayes  rule 
E(F(t)  I  kj . it,  ) .  The  former  evaluation  requires  only  partial  specification 


This  approach  of  approximating  Bayes  rules  in  some  nonparamecric  problems 
by  Bayes  estimators  within  the  class  of  linear  combinations  of  a  given  set  of 
functions  on  the  sample  space  has  been  proposed  by  Goldstein  (1975).  He 
applies  his  result  to  derive  linear  polynomial  estimators  of  the  mean  and 
moments  of  the  unknown  distribution  function.  However,  Goldstein  does  not 
specifically  treat  the  bioassay  problem.  The  objective  and  loss  function 
studied  here  are  also  different  from  Goldstein’s.  Integrated  loss  function  is 
used  here,  since  the  objective  is  to  estimate  the  whole  distribution  function. 

To  use  Ferguson's  process  indexed  by  a,  the  statistician  needs  to  specify 


c .  The  measure  m(-=>,t]  can  be  rewritten  as  MF^(t),  where  M  =  sC-®,00)  and 
F  \t)  =  a  ( -* ,  t  j ,  >1 .  It  is  known  that  EF(t)  =  Fn(t)  and  var(F(t))  = 

W  U 

F  . t ;  ( 1-F  ( t) )  /  (M-r  1 )  where  F  is  the  unknown  random  distribution  chosen  by  the 


Ferguson's  process.  The  function  F.  represents  the  statistician's  prior  belie 


'h  the  shaps  of  F,  and  V.  represents  the  degree  of  concentration  of  F  around 


ates  that  F  concentrates  more  around  F,. 


4 


The  linear  Bayes  rule  d  derived  in  this  paper  can  be  summarized  as 

follows.  For  t  to  be  one  of  the  stimulus  levels,  say  t  «  t.,  the  linear 

J 

A 

Bayes  rule  d(t^.)  is  a  linear  combination  of  Fq(z^),  and  k^/n^  -  F^Ct^)  for 
i=l,...,L.  The  coefficients  are  ratios  of  two  determinants  of  LxL  matrices 
where  the  entries  are  functions  of  M,  Fq,  and  n^,...n^.  For  values  of  t  not 
at  the  drug  levels,  say  for  example  t^  <  t  <  d(t)  can  be  interpolated 

from  d(t^)  and  d(t  using  F^.  Perhaps  the  most  interesting  result  of  this 
paper  is  this  interpolating  formula.  It  reduced  the  problem  of  evaluating 

A  A  A 

d(t)  for  each  t  to  the  problem  of  just  evaluating  L  of  them,  d( t j ) , . . . ,d( t  ) . 

One  of  the  drawbacks  of  the  linear  Bayes  estimator  is  that  it  may 
not  be  nondecreasing.  Fortunately,  nonmonotonicitv  does  not  occur  very  often 
If  monotonic  restriction  is  a  great  concern  to  the  user,  the  pool-adjacent- 
violators  algorithm  (Barlow,  Bartholomew,  Bremner  and  Brunk,  1972,  pp.  13-18) 

A  A 

car.  be  used  for  obtaining  the  isotonic  regression  on  d( t  ) ,  . . .  ,d( t  ) . 

1  L 

Nevertheless,  this  rule  is  not  Bayes  to  the  original  problem.  It  is  also  not 
clear  what  to  do  for  t  not  at  the  drug  levels.  An  alternative  formulation  of 
linear  Bayes  estimation  under  order  restriction  might  be  of  some  interest  and 
could  be  explored  later.  The  present  paper  adopts  the  one  without  an  order 
restriction  for  computational  simplicity.  This  formulation  enables  us  to 
cerive  the  interpolating  formula  given  in  Theorem  2  for  the  linear  Bayes  rule 
The  general  discussion  on  the  Dirichlet  process  and  the  mathematical 
model  for  the  quantal  bioassay  is  given  in  Section  2.  The  derivation  of  the 
linear  Bayes  estimator  is  given  in  Section  3.  Some  asymptotic  results  and 

ii..iiU.V5  uc..3viui  Zu  C  Vv*  vj  • j  C  I  u  C  l  lie  OdVeS  rules  uTc  i.fi 

Sections  4  and  5.  A  numerical  example  comparing  the  Bayes  estimators  and 
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linear  Bayes  estimators  for  L  -  2  is  given  in  Section  6.  A  SAS  program  is 
given  in  the  appendix  to  illustrate  the  ease  of  computing  linear  Bayes 
estimators  in  practice.  The  program  computes  the  linear  Bayes  estimators  for 


all  possible  observations  for  n^  *  n^  =  3.  The  PROC  MATRIX  procedure  which 


has  the  capability  of  inverting  a  matrix  is  used.  Another  numerical  example 
given  in  Section  7  using  the  data  in  Cox  (1970,  p.9)  illustrates  the  use  of 
linear  3ayes  estimators  in  practice. 


2.  PRELIMINARY  WORK 

A  brief  description  of  the  Dirichlet  process  and  the  quanta!  bioassay 
problem  is  given  here. 

Ferguson  (1973)  constructs  a  prior  on  the  space  of  distribution 
functions.  The  prior  indexed  by  a  chooses  a  random  distribution  function  ?. 
This  is  denoted  by  ?-:  D  (a) .  The  prior  can  be  defined  as  follows.  Given  any 
measurable  partition  3  ....  ,3,  of  R1 ,  the  joint  distribution  of  the  random 

1  it 


probabilities  ?( 3. ) , . . . ,P( 3.  )  is  a  Dirichlet  distribution  with  parameters 

1  K. 


csCSj)  , . . .  .2(5^) . 


Giver,  the  loss  function 


L(F,F)  =  /  (F(c)  -  F( t ) ”  dW(t), 


it  is  shown  by  Ferguson  that  the  Bayes  estimator  of  F  for  the  no  sample 

*.  O  *  -  C  —  x  S 


•(t)  =  EF( t )  =  o(-»,tj/o(R)  =  MrQ(t)/M  =  FQ( t)  . 


The  Haves  estimator  of  F,  when  we  observe  a  sample  x,,...,x  from  F,  is 

i  n 


given  oy 


6 


f; 


I 


Fn(t)  =  E(F(t) |x^ . xr)  =  {a(-®,t]  +  E  I(xi<t)}/(M+n) 

1  k 


TT-  FA(t)  +  r~  F  (t)  , 
M+n  0  M+n  n  ’ 


where  F^  is  the  sample  empirical  distribution  function. 


Recall  the  bioassay  problem,  is  the  number  of  subjects  react 
positively  at  the  stimulus  level  t^.  Therefore,  the  likelihood  function  is 
given  by 

L 


L(F)  =  n  ("M  {F(t.)}Ki  {  1  — F  ( t .  )}ni-Ki  . 
i-1  Ki'  1  L 


k-f 


ni"ki 


(2.1) 


Given  the  Ferguson's  prior  F-.D(a)  and  L  =  2,  it  is  shown  by  Antoniak  that  the 
posterior  distribution  of  F  is  a  mixture  of  Dirichlet  processes.  Define 
“1  =  =  -■KF(j(t,)-F0(ti)),  and  »  M(l-Fg(t?)).  It  is  shown  by 


Antoniak  that  the  posterior  distribution  of  FCt^),  Fft^-FCt^),  is  given  by 


-  *  *  "  ^ 


i=:  :=o 


aij  D  (  =  1+ki'i'i>  P 7+ni-ki+k2~i_^  *  e3+n2_k2+'’ )  ’ 


(2.2) 


where  the  mixing  indices  are 


k2  nj-kj 

a..,  =  b .  /Z  T.  b. 


1J  i*0  >0 


and 


n  f  -k  1  w k2 -v  r(gi+k;-^i)  r(3?+n|-ki+k?-i--i)  f  ( g ^+n?-k?-rj  ) 


1 ."i)f  2) 

■  j  J'l  1 


?(sl)  r(g2)  r(?3) 


Antoniak  has  also  given  a  good  heuristic  explanation  for  the  mixture.  While  it 
is  known  mat  <,  ocservacions  fall  in  (0,c_j  ana  n  -k.  observations  tali  in 

~  4m  11 

(tpl],  we  don't  know  how  many  of  the  k^  fall  in  ( 0 ,  t  ^  ]  and  how  many  of  the 

n  -k  fall  in  (t,,lj.  Let  i  denote  the  number  of  observations  in  (0,t,] 
i  1  <.  1 

arising  from  k7,  and  let  j  denote  the  number  of  observations  in  (t^.l)  arising 


WV| 


from  n^-k^ 


Since  i  and  j  are  unknown,  the  posterior  distribution  is  a 


mixture  over  all  possible  values  of  i  and  j 
the  expression  (2.2). 

The  Bayes  estimators  can  be  obtained  as 


FCtj) 

k2 

=  Z 

n  l  ~k  i 

Z 

a.  . 

i-0 

L_i. 

II 

O 

ij 

F(t,) 

k2 

=  Z 

nl~kl 

z 

a .  . 

2 

i=0 

3=0 

ij 

given  by  the  double  summation  in 


3)  +  ki  + 


M 


ni  +  n2 


Si  +  Bo+ni-rko-j 
M  t  nj  +  no 


When  L  =  3  ,  it  takes  six  summations  to  describe  this  mixture.  In  genera 
it  cakes  L(L-l)  summations  to  describe  the  mixture.  To  circumvent  the 
difficulty  in  computing  the  Bayes  estimators,  the  linear  Bayes  estimator  is 
proposed  and  derived  in  the  next  section. 


3.  DERIVATION  OF  THE  LINEAR  BAYES  ESTIMATOR 

2 

liven  the  loss  function  L(F,d)  =  /(F(t)-d(t))  dW(t),  the  linear  Bayes 

o 

estimator  is  obtained  by  minimizing  the  function  E( F( t )-d( t ) ) “  for  each  t. 

Let  us  first  restrict  our  attention  to  estimating  F(t)  at  the  stimulus  levels, 

namely  when  t  *  t^,  j=l,...,L.  The  following  Theorem  1  says  that  d ( c ^ )  is 

actually  a  sum  of  F„(c.)  and  a  linear  combination  of  k./n.-  F„(t.),  i=i,...,L, 

0  j  l  l  0  l 

where  the  coefficients  defending  on  t,,...,tT  can  be  computed  from  the 

1  L, 

determinants  of  certain  covariance  matrices. 

We  need  the  following  two  notations  for  Theorem  1.  Let  cov(k)  denote  the 

covariance  matrix  of  k^,...k  and  let  A(i,t)  denote  this  matrix  with  the  ic^ 

T 

column  replaced  by  lcov(k^  ,F(t) ) ,  cov(k,  ,F' t) ) , . . .  ,cov(  rC^,F(.t) )  ]  ,  namely 


, -c. 


'*1 


var  k, 

i 


cov( k, , k_.) 
x  2 


cov(k^,k9)  var  k„ 


1  column 
cov( k ^ ,F(t)) 

cov(k2>F(t) ) 


cov( k j , k  ) 


cov( k^ ,k  ) 


cov(kL>F(t))  ...  var  k^ 


>_reu_l_:  In  the  bioassay  problem,  let  us  assume 

(1)  the  potency  curve  F(t)  is  a  random  distribution  chosen  from  the 

Lchiet  process  with  parameter  MF.(t)  (-®<t<®); 

L) 

(2)  n,,...,n.  subjects  are  created  at  stimulus  levels  t  <t0<...<t 
lettively,  and  k,,...,k  are  the  observed  numbers  of  subjects  that  react 
Ltively  at  each  level; 

(3)  the  loss  function  is  L(F,d)  =  J(F( t)-d( t) )“  dW(t). 


.  w  WU1  uCu-bi'VU 


ru*65  co  cne  cj.&ss  ox  i.  i  n  ti  3.  r  conioinacions  or 


, k.  ar.c  1  for  each  t,  then  the  Bayes  rule  in  this  class,  evaluated  at  the 
.us  levels  t.,j=l,...,L  is  given  as  follows: 


d(t  )  =  F  (t  )  +  l  n.k.(j)  (k./n.  -  F  (t  ; 
J  J  1=  1 


Xi(j)  =  |A(i,tj ) j/|cov  k| 


j* I  denotes  the  determinant  of  matrix 


(3.1) 


(3.2) 


-n  »» .  * 


PROOF .  Given  t  =  t.  ,  we  seek.  Xq( j ) ,  • • • ,X^( j)  to  minimize 


b(X)  =  ECF(c  )  -  X0(j)  -  z  X1(j)  k.) 


1=  i 


where  E  denotes  the  expectation  over  the  joint  distribution  of  the  observations 
and  the  parameter  F(t,).  Differentiating  <$(\)  with  respect  to  Xq( j ) , . . .  .X^C j  ) 
respectively,  we  obtain  L+l  normal  equations: 

T 

„  u 

E(F(tj)  -  XQ(j)  -l  X.(j)k.)  k^  =  0,  1-0,..., L, 

where  k.  =  1. 


Therefore,  from  the  equation  for  1=0,  we  obtain 


X^(j>  =  EF( t . )  -  )  X  .  ( j )  E(k.) 

i-1  1 


(3. 


Substitut 
and  solvir 


ing  this  expression  for  X..  into  each  of  the  L  equations  with  1  =  1 . L 


:g  the  system  of  L  linear  equations,  we  obtain 


X.(j)  =  i A(i,t  ) j/ jcov  k(  i=l,...,L. 


M  r* 


?/sr,  f rod  3  c  o  3  -  i  o  r.  (3*3^ 


X0(j)  -  ^•iO)niFu(ti), 


d(t J  =  XQ(j)  +  j  Xi(j)ki 
1  =  1 


=  F.(t.)  * 


n.X.(j)  (k./n.  -  F ,  ( t . ) ) . 


Next,  we  want  co  estimate  F(t)  for  t  other  than  the  drug  levels,  for 


example  t.<t<t.  We  can  find  d(t)  for  each  t  by  the  same  method  as  in 
J  J+1 

Theorem  1,  where  t^  is  replaced  by  t  in  equations  (3.1)  and  (3.2).  If  d(t)  can 
easily  be  derived  from  Fq,  d(t^),  and  d(t^.  +  ^),  then  time  and  effort  can  be 
saved  in  estimating  F(t),  since  we  no  longer  need  to  compute  the  determinant  of 
A(i,t)  for  each  t.  The  following  theorem  enables  us  to  linearly  interpolating 
d(t)  by  using  d(tj)  and  d(t^+^). 

Theorem  2:  With  the  assumptions  of  Theorem  1,  and  restriction  of  decision  rules 
to  linear  combinations  of  k  ,...,k  ,  and  1,  the  Bayes  rule  for  F(t)  is  given  by: 

F.(t.  ,)-F.(t)  „  F.(t )  -  F.(t . ) 

d(t)  =  ^7^ —  d(t4)  f  yt - ~jr)d(t..+1) 

‘ 0  w j+l  F0vtjy  J  F0(tj+1)  *ouj)  J  +  i 

for  t.  <  c  <  t  j,  j=0,...,L  with  FQ(tQ)  -  0  «  d(tQ), 

F0^CL-i-P  =  1  *  a(CL+l^' 


PROOF .  By  the  same  method  as  Theorem  1,  we  obtain 


d(t)  =  FQ(t)  *  J  ni  Xt(t)  (ki/ni  -  FQ(ti)) 


-  X0(t)  +  E  \.(t)k., 
i=  1 


(3.6) 


where  \,(t)  =  F  (t)  -  £  \.(t)  n_^  F^t^,  and  \A  t)  ,  i*l , . . .  ,L  is  obtained 
from  the  right  side  of  the  equation  (3.«)  with  F(t^)  replaced  by  F(t).  Note 
that  the  entries  cov(k_.,F(t))  can  be  computed  similarly  as  in  equation  (3.5). 
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cov(k.,F(t))  -  :  n.  F0(ti)(l-F(t))/(«+l),  if  t.  <  t, 
n.  F0(c)(l-F0(c.))/(M+1),  if  tL  >  t. 

For  t.<  t  <  t  j,  it  can  be  verified  chat  for  i*l,...,L 


,/o(ci+i)  f0(C) 

,--Cc))=  cov(s  ,  ■  ( c  .  )h : - ^ - 

i.  1 


F0(t)'F0(Ci) 

+  cov(kifF(cj+1))  —-——1— 


‘o'  j+r  ‘O'  j- 


O'  j+l7  0  j 


.  (3.7) 


Then  by  a  straightforward  matrix  multiplication  using  equations  (3.4),  (3.6) 
and  (3.7),  we  can  show 


F,(ti_1)-F0(c)  ^ 

d<t:)  *  tjt  “'U-i' 


F„(t)-F„(t.) 


_ 3. 


d(c  .  ,). 


j  ■  j-r  j  j 


4.  ASYMPTOTIC  RESULTS  FOR  d(c). 

*  A 

ay  evaluating  ^(j),  we  are  able  to  show  d(t  )  is  an  asymptotically  unbiased 
:  consistent  estimator  of  F(c^). 

^rer  3:  Given  the  assumptions  of  Theorem  1,  let  ®  for  all  i*l,...,L.  The: 


t4).ri  •*  ■(;.),  and  d(t.)  is  a  weakly  consistent  estimator  of  F(t.)  for 
-  j  J  J 


■Or .  From  equation  (3.1),  we  have: 


d( 


t  )  ■=  F.^(t  )  +  T  n  A  (j)(k  /n  -  F  (t  )). 
3  u  3  i-111  1  u  1 


(4.1) 


It  can  be  verified,  by  evaluating  equation  (3.2),  that 


0(l/n  ) 


if  i*j. 


n.X.(j)  -  J 


1  +  T  0( I/n  )  if  i  *  j . 
i-1 


Tneref ore. 


E(d(t.)|F)  -  F0(t.)  +  l  n.Xi(j)(F(t.)  -  FQ(t.)) 

L 

“  FQ(t  )  +  (1  +  l  0(l/n  )|(F(t  )  -  F  (t  )) 
j  X  J  UJ 


+  y  0( l/n  ) ( F C t  )  -  F  (t  )) 
i*j  J 


F(t..)  as  a.  +  «  for  i»l,...,L. 


Then,  from  equations  (4.1)  and  (4.2), 


ft  .  1 

var(d(  t  )  ;  F)  =»  E  £  n  X  (j)((k  ,'n  )  -  F(t.)) 

J  i=l 


I  (j)4  var(k  /n  ) 

i=l  1  1 


II  n«MJ>  n.X.Cj)  cov(k./n. ,  k/n) 
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JOT 

wv 


L 

Z 

i-1 


(1  +  Z  0(l/ni))^  .  F(tj)(l-F(e  ))/n 


$ 


+  l  0( l/n. ")  .  F(c.)(l  -  F( t  ) )/n 
,  .  J  i  i  i 

i*j  J 


+  0  as  n.  +  «  for  i=l . L. 

i 


(Note  that  given  F,  is  independent  of  k^.) 


Therefore, 


d(t^)  is  a  consistent  estinator  of  F(t^). 


5.  LIMITING  BEHAVIOUR 


(5.1) 
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This  is  verified  in  a  similar  manner  to  (1). 

(3)  If  i“ j ,  then  Xj(j)  =  j A( j , t^ ) | / j cov  k|. 

When  M  -*■  0,  we  can  factor  n,  from  the  column  of  the  denominator,  and  show  the 

J 

A 

ratio  of  the  two  remaining  determinants  -*•  1.  Therefore,  lim  \.(j)  =  1/n.. 

tt*0  J  J 

From  equation  (3.1),  we  have 


lim  d(t.)  =  lim  [F  (t.)  +  £  n.\.(j)(k 

.  M->0  J  M+0  J  i=l 


/n.  -  F0(t.)J 


F0(t.)  +  (k./n,  -  F0(tj))  -  k./n.. 


REMARK: 


Antoniak  has 


given  an  example  for  M  ♦  0  with  n^  =  n?  =  100,  k^  =  1 ,  k^  =  99, 


t,  =  1/3,  t0  =  2/3,  and  F  (t)  =  t  ■=  £0,1]  .  He  obtains  the  Bayes  estimate 

i  i.  u 

F'i/3)  =  F( 2/3)  =  1/2  and  F(t)  is  piecewise  linear  on  10,1]  (F(0)  =  0,  F(l)  -  1). 

A  A 

The  linear  3aves  estimate  in  this  case  is  given  bv  df  1/3)  =  1/100,  d(2/3)  = 

*  •  \  y  7  v.  j 

A  A  A 

99/100  and  d(t)  is  piecewise  linear  (d(0)  =  0,  d ( 1 )  =  1).  Paradoxically,  the 
linear  Bayes  estimate  d  in  this  case  is  intuitively  more  reasonable  than  the 
usual  3aves  rule  F. 


6.  COMPARISON  OF  LINEAR  BAYES  TO  BAYES  ESTIMATES 
Comparisons  between  Bayes  and  linear  Bayes  estimates  are  made  in  this 
section.  For  L=2,  the  Bayes  estimator  of  F  can  be  evaluated  as  in  Section  2. 
Let  us  assume  n^=  n^=  3,  i.e.  three  subjects  are  tested  at  each  stimulus  level. 
The  linear  Bayes  estimates  and  the  Bayes  estimates  evaluated  at  the  two  drug 
fcvcfs  with  ditiereut  M  are  given  m  Table  1  and  Table  2  iur  an  the  possible 
values  of  k^and  k^.  The  prior  probabilities  for  each  of  the  events  are  also 


;***> 


/i 
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given  in  the  last  column.  The  prior  probabilities  are  evaluated  by  E(L(F)) 
with  L(F)  given  by  equation  (2.1)  with  appropriate  n, ,  n_,  k.  and  k_.  Both 
cables  use  the  same  Fq,  standard  normal  distribution,  M= 1  and  M*10  are  used 
respectively  in  Table  1  and  Table  2.  The  stimulus  levels  are  administered  at  -1 
and  1.  It  can  be  seen  that  the  linear  Bayes  estimates  approximate  the  Bayes 
estimates  fairly  well  for  M=10.  For  M=l,  the  approximation  is  satisfactory  for 
cne  cases  k,  <  k., .  The  worst  case  occurs  for  the  case  k,  =»  3  and  k,  »  0. 
Fortunately,  che  prior  probability  for  this  event  to  occur  is  only  0.0002. 

7.  NUMERICAL  EXAMPLE 

The  data  for  this  example  are  taken  from  Cox  (1970  p.  9).  In  all,  150 
subjects  are  tested  at  5  different  concentrations  of  the  stimulus  with  30 
subjects  at  each  of  the  dosages.  The  numbers  of  deaths  are  recorded 
respectively  at  the  dosage  levels.  Three  different  priors,  X(2,0.5),  X(2,l) 
and  with  various  M  are  used  for  computing  the  linear  Bayes  estimates 

evaluated  at  tne  rive  dosages. 

It  can  be  observed  chat  the  linear  Bayes  estimators  are  not  monotonic  in 
M.  As  M  increases,  the  linear  Bayes  estimates  approach  the  prior  guess 
F  23  expected.  When  M=l,  the  linear  Bayes  estimates  are  relatively  insensitive 
to  the  choice  of  che  prior  guess.  They  are  all  close  to  the  isotonic  maximum 

x  1  II  TiCCC 


8.  CONCLUSIONS 


The  nonparametric  Bayesian  method  is  applied  to  estimate  the  potency 
curve  in  the  quantal  bioassay  problem.  To  circumvent  the  computational 
difficulty  in  evaluating  the  Bayes  rule,  the  linear  Bayes  estimator  is 
proposed.  The  linear  Bayes  estimate  can  be  obtained  easily  by  linearly 
interpolating  the  estimates  evaluated  at  the  dosage  levels.  To  obtain  the 
estimates  at  the  dosage  levels,  it  is  required  to  invert  a  LxL  matrix  which 
can  be  done  by  many  standard  statistical  packages.  Some  numerial  examples  are 
given  to  demonstrate  the  ease  of  obtaining  the  linear  Bayes  estimate. 
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Table 

1.  Comparison 

between 

linear  Bayes 

and  Bayes 

estimates 

h-i  r0(t 

!>  -  °‘: 

1587  F0(t2) 

-  0.8413 

Prior 

Observations 

Linear  Bayes 

Bayes 

Probability 

i. 

k2 

d(c1 ) 

d(t2) 

FCtj) 

F(t,) 

pr 

0 

0 

0.0101 

0.2089 

0.0227 

0.1964 

0.0599 

0 

l 

0.0221 

0.4572 

0.0372 

0.4422 

0.0632 

sj 

0.0342 

0.7055 

0.0419 

0.6975 

0.0925 

J 

3 

0.0462 

0.9538 

0.0415 

0.9585 

0.5180 

1 

>j 

0.2584 

0.2209 

0.1655 

0.3138 

0.0049 

1 

1 

0.2704 

0.4692 

0.23oo 

0.5029 

O.uC'92 

l 

2 

0.2325 

0.7175 

0.2855 

0.7145 

0.0163 

i 

3 

U.2945 

0.9658 

0.3022 

0.9581 

0.0925 

- 

0 

0.5067 

0.2330 

0. 3084 

0.4313 

0.0012 

1 

4. 

1 

0.5187 

0.4813 

0.4085 

0.5916 

0 . 0  o  o  i 

/ 

-> 

0.5303 

0.7296 

0.4971 

0. 7 632 

0.0092 

2 

3 

0.5428 

0.9779 

0.5578 

0.9b23 

0.0632 

3 

0 

0.7  550 

0.2450 

0.4512 

0.5^83 

0 .  *  j 1  J  C  2 

j 

1 

0.7670 

0.4933 

0. 5687 

0.6916 

0.0012 

3 

n 

4m 

0.7791 

0.7416 

0.6862 

0.8345 

0.0049 

3 

3 

0.7911 

0.9899 

0.8036 

0.9773 

0.0599 

Table  2. 

Comparison  between 

linear  Bayes 

and  Bayes 

estimates 

) 

) 

) 

1 

M-10  F 

t 

0(tl)  “  °* 

1587  Fq(c2 

)  =  0.8413 

* 

1 

Prior 

j  Observations 

Linear 

Bayes 

Bayes 

Probability 

\  ki 

\ 

k2 

d(c1) 

d(t9) 

FUj) 

F(t2) 

Pr 

|  u 

» 

0 

0.0941 

0.6431 

0.0992 

0.6380 

0.0078 

1 

1  0 

1 

0.1052 

0.7195 

0.1087 

0.7161 

0.0523 

- 

0. 1 lb4 

0. 7 960 

0.1178 

0.7946 

0.1824 

0 

3 

0.127b 

0.872** 

0.1264 

0.8736 

0.3823 

* 

•j 

0. 1705 

0.6543 

0.1617 

0.6631 

0.0028 

i 

i 

1 

0.1817 

0.7307 

0.1767 

0.7357 

0.0207 

1 

* 

') 

0.1929 

0.8071 

0.1913 

0.8087 

0.0797 

l 

3 

0.204Q 

0.8836 

0.2054 

0.8822 

0.1824 

■) 

i 

*> 

0.2*69 

0.6655 

0.2242 

0.6882 

0.0006 

\ 

1  : 

■ 

i 

0.2581 

0.7419 

0.2444 

0.7556 

0.0048 

* 

r 

t 

0 

0. 2b93 

0.8183 

0.2643 

0.8233 

0.0207 

i 

3 

0.2805 

0.8943 

0.2839 

0.8913 

0.0523 

;  3 

a 

0 

0.32  3** 

0.6766 

0.2867 

0.7133 

0.0001 

• 

3 

» 

1 

0.3345 

0.7531 

0.3118 

0.7758 

0.0006 

> 

k  .  3 

1 

0.3457 

0.8295 

0.3369 

0.8383 

0.0028 

!  3 

3 

0.3569 

0.9059 

0.3620 

0.9008 

0.0078 

Table  3.  Linear  Bayes  estimates  for  Cox’s  data 


Log^vConcentrat 

ion)  0 

1 

2 

3 

4 

of  deach 

2 

8 

15 

23 

27 

n 

i 

0.067 

0. 2o7 

0.500 

0.767 

0.900 

~  *  m  5\ 

V\  As  — 

V 

0.0000 

0.0228 

0.5 

0.9773 

1.000 

a 

A 

A 

A 

_ 

M 

d  ( 0 ) 

d  ( 1 ) 

d  ( 2 ) 

d  ( 3 ) 

d  ( 4 ) 

l 

0.0645 

0.2606 

0.5005 

0.7706 

0.9032 

0.0334 

0.1594 

0.5022 

0.8491 

0.9499 

1  00 

0 . 0 1 5  5 

0.0894 

0.5U07 

0.9125 

0.9769 

3 -‘J 

O.OUoi 

0.0497 

0.5002 

0.9506 

0.9909 

5  2  0 

0.0038 

0.0397 

0.5001 

0.9605 

0.9943 

r  s  v  f  i  in 

*  -  A  s  —  *  ^  y 

t,) 

o  ,  \  i  o  a 

v  •  v  ^  —  O 

0.1587 

0.5 

0.8413 

0.9772  ! 

i 

0 . 0636 

0.2633 

0.5007 

0.7676 

0.902a 

J'j 

0.0476 

0.2149 

0.5026 

0.7939 

0.9375 

A  V  J 

0.0349 

0.1854 

0.5011 

0.8172 

0.9584 

3  3  3 

0.0277 

0.1694 

0.5003 

0.8313 

0.9697 

—  > 

0.0259 

0. 1  65-+ 

0. 5uo2 

0.8350 

0.9725 

:  i; 

0.3085 

0.4013 

0.5 

0.5987 

0.69151 

J 

i 

0.0794 

0.2692 

0.5019 

0.7576 

0.S922 

J  ■  J 

0.2115 

0.3380 

0.5020 

0.6667 

0.7329 

100 

0.2675 

0.3741 

0.5008 

0.6276 

0.7307 

300 

0.2931 

0.3910 

0.5003 

0.6097 

0*7  Qo4 

5o0 

0.2990 

0.3949 

0.5002 

0.6055 

0.7007 
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APPENDIX:  SAS  source  program  for  computing  the  linear  Bayes  estimates, 

DATA  OBS; 

INPUT  01-02; 

CARDS; 

0  0 
0  1 
0  2 


DATA  SIZE; 
INPUT  S1-S2; 
CARDS; 

3  3 


DATA  PARAM; 

Pl-PROBNORM(-l); 

P2=?R0BN0RM(  1 ) ; 

M-l; 

PROC  MATRIX; 

FETCH  0  DATA-03S; 

FETCH  N  DATA-SIZE; 

FETCH  B  DATA-? ARAM; 

M=B(  1,3); 

T  —  ?  . 

—  » 

DO  K»1  TO  16; 

X-J(L,L); 

C=L(L ,L) ; 

A-J(L,L); 

LAMDA-J ( 1 ,L) ; 

EST*J( 1 ,L) ; 

DO  1=1  TO  L; 

X(I,I)=N(K,I)*(M+N(K,I))*3(1 ,1)*(1-B(1 ,1))v‘/(M+l); 
END , 


a 


APPENDIX:  Contiaed . . . . 


k 


E 


DO  1*2  TO  L; 

DO  J=1  TO  I- I; 

X(I,J)«N(K,I)*N(K,J)*B(I ,I))#/(M+1); 
X(J,I)=X(I,J); 

END; 

END; 

DO  J*I  TO  L; 

DO  1=1  TO  L; 

IF  I<=J  THEN  DO; 

C(I,J)=N(K,I)*B(1 ,1)*(1-B(  1  ,J)  )•!/  (M+l) ; 

END; 

ELSE  DO; 

C(I,J)=N(K,I)*B(1 ,J)*(1-B(1 ,I))#/(M+1); 

END; 

END; 

END; 

DO  J=1  TO  L; 

A(  ,  J)=INV(X)*C( ,J) ; 

3B=3( 1 , l  2); 

LAMDA(  1  ,  J )  =B(  1 ,  J  ) -( AC  ,  J  )  ,v‘N ( K ,  V  )  '  *BB  ' ; 

ESTCl ,J)=A( ,J) '*0(K, ) ’+LAMDAC 1 ,J) ; 

END; 


